TELKOMNIKA Telecommunication Computing Electronics and Control 
Vol. 21, No. 1, February 2023, pp. 123~132 
ISSN: 1693-6930, DOI: 10.12928/TELKOMNIKA.v21i1.20746 o 123 


Hippocampus’s volume calculation on coronal slice’s for 
strengthening the diagnosis of Alzheimer’s 


Retno Supriyanti, Yogi Ramadhani, Eko Wahyudi 


Electrical Engineering Department, Faculty of Engineering, Jenderal Soedirman University, Purbalingga, Indonesia 


Article Info 


ABSTRACT 


Article history: 


Received May 25, 2021 
Revised Nov 04, 2022 
Accepted Nov 14, 2022 


Alzheimer’s is one of the most common types of dementia in the world. 
Although not a contagious disease, this disease has many impacts, especially 
in socio-economic life. In diagnosing Alzheimer’s and using interview 
techniques, physical examination methods are also used, namely using an 
magnetic resonance imaging (MRI) machine to get a clear image of the 


patient’s brain condition, with a focus on the hippocampus and ventricular 
area. In this paper, we discuss the calculation of the volume of the 
Keywords: hippocampus, especially the coronal slice, to provide information to doctors 
in making decisions on diagnosing the severity of Alzheimer’s. Using the 
basis of volume calculations, we made a 3D visualization reconstruction of 
the coronal hippocampus slice area in order to make it easier for doctors to 


Alzheimer’s diagnosing 
Coronal slice 


Hippocampus analyze the condition of the hippocampus area, which in the end will be used 
MRI as a recommendation in the classification of the severity of Alzheimer’s. 
Severity Our experimental results show, the lower the severity, the bigger the volume, 


the more slices, and the longer the counting time. 
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1. INTRODUCTION 

Alzheimer’s disease has a record of an increase in sufferers every year; this disease will 
undoubtedly attack most people who have the elderly. However, this fact is inversely proportional to the 
availability of services provided in several health agencies. It is considered that this disease does not 
provide high profits, such as cancer, heart disease, kidney disease, and stroke, which are generally suffered 
by most people in developing countries like Indonesia. Many people with Alzheimer’s are ignored. This 
disease will cause the sufferer to experience memory impairment, impaired judgment, decision-making, 
overall physical orientation, and problems with speech. This disease is divided into two subtypes based on 
age, namely early-onset Alzheimer’s disease (EOAD) and late-onset Alzheimer’s disease (LOAD). EOAD 
accounts for about 1% to 6% of total cases and ranges from 30 years to 60 or 65. However, LOAD is the 
most common form of Alzheimer’s affecting humans aged 60 and 65 years. Both EOAD and LOAD are most 
likely to occur in people with a family history of Alzheimer’s disease. As many as 60% of EOAD cases are 
from families with a history of Alzheimer’s disease [1]. Several stages are carried out in the medical field to 
diagnose Alzheimer’s, namely history, clinical examination, and magnetic resonance imaging (MRI) 
radiological support. Diagnosis based on the patient’s behavior. 

In radiology, MRI examination of the observed part is the hippocampus. Alzheimer’s disease will 
experience abnormal conditions in the hippocampus. Neurodegeneration occurs in the brain, and oxidative 
stress is associated with an increase in amyloid beta (Af) deposits. AB causes H202 accumulation in 
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hippocampal neuron cultures and neuroblastoma cultures. There is a shrinkage in the hippocampus according 
to the increase in the critical dementia ratio (CDR) value in people with Alzheimer’s [2]. In radiological 
examinations, the most common observation of MRI images by radiologists is by photo printing. For simple 
cases, at least two photo printing sheets must be done on each axial, coronal, and sagital cut. In complex 
cases, there may be more than five photo prints on each piece for observation. In low-resolution MRI 
machines such as those available in health services in developing countries such as Indonesia, this MRI machine 
has a weakness. The frequency of use is limited because machine precision needs to be maintained to decrease 
performance. The processing process to get the desired organ image takes a long time. For example, in the most 
straightforward cases, such as slight scratching of the bone, it takes up to half an hour, and in complex cases, 
it takes more than an hour. Every time the contrast improvement in the image is made, an additional cost must 
be added. This long-term research aims to develop a simple system that can be used to help doctors diagnose 
Alzheimer’s patients based on severity. We have done some preliminary research before [3]-[6] but the 
research we did still covers the area that has not yet reached the volume calculation. 

There is much research on Alzheimer’s and its relationship to image processing. Among them are 
the following. Keserwani et al. [7] classified Alzheimer’s based on the texture analysis of the hippocampus 
area. The hippocampus area was extracted using Gabor segmentation. Texture features are represented on a 
slice-by-slice basis with mean and the standard deviation of the magnitude of Gabor’s response. Ozsahin et al. [8] 
uses backpropagation neural networks for automatic classification of 4 patient groups, including Alzheimer’s 
disease, late mild cognitive impairment, early mild cognitive impairment, and significant memory concern, versus 
standard control on 18F-florbetapir positron emission tomography (PET) images from the Alzheimer’s disease 
neuroimaging initiative (ADNI) database, with 100 images for each group. Bothraa et al. [9] developed an 
automated technique to detect Alzheimer’s symptoms using an experimental approach. They divided their 
research into two parts, first is to get a dataset of ADNI disease and perform image processing on it that will 
be used to train the system. Furthermore, using the deep learning algorithm to detect a disease from this 
neuroimaging data. Giang et al. [10] they proposed using a fast-multiple kernel learning framework, known 
as fast multi kernel learning framework-dimensional reduction ((MKL-DR), to classify Alzheimer’s disease 
(AD) patients and cancer patients. Maqsood et al. [11] developed an automated dementia detection and 
classification system. The system classifies dementia patients and identifies the four stages of dementia 
development. they utilize transfer learning to classify images by refining pre-trained convolutional network, 
AlexNet. Pelsmacker et al. [12] explored the influence of message as a mediator in the processing of 
Alzheimer’s consciousness with various arguments, strength, and valence. Islam et al. [13] developed a 
classification method for 3D organ images are rotation and invariant translation. They extract a representative 
two-dimensional (2D) piece along the plane best symmetry of 3D images. furthermore, they use these slices 
to represent 3D images and uses a 20-layer deep convolutional neural network (DCNN) to perform 
classification tasks. Bhagyashree et al. [14] researched by developing software to link traumatic brain injury 
to the severity of Alzheimer’s. Park and Moon [15] reviewed the use of MRI as a basis for supporting image 
processing-based Alzheimer’s diagnosis. Cuthbert et al. [16] classified the severity of dementia using transfer 
learning with visual geometry group model 16 (VGG16) using fast AI. Kleij et al. [17] conducted a study on 
the accuracy of measuring brain parenchymal volume (BPV) and intracranial volume (ICV) by performing 
three-dimensional segmentation. Hadi et al. [18] proposed region-based segmentation of brain tumors using 
an active contour model. Wu et al. [19] proposed a three-dimensional network model for brain tumor 
segmentation by adopting an encoder-decoder structure and replacing convolution with group convolution to 
improve network performance. This paper will discuss calculating the hippocampus’s brain volume based on 
the coronal slice of MRI segmentation images based on time, the value of CDR, gender, and age. 


2. RESEARCH METHOD 
2.1. Input data 

In this research, we used input data from the open access series of imaging studies (OASIS) 
database. The OASIS project aims to create a set of MRI data sets that are freely accessible and free for the 
scientific community [20]-[21]. Figure 1 shows a sample of a coronal slice of hippocampus image data from 
Alzheimer’s patients based on the CDR 0 to 2. A value of 0 indicates no dementia, 1 is mild dementia, and 2 
is moderate dementia. 


2.2. Active contour segmentation 

This paper is a continuation of our previous research [4]-[6], [9], [10], where the segmentation used is the 
active contour segmentation discussed in the previous paper. Active contour is a segmentation method using a 
closed curve model that can move broader or narrower. Active contour can move broader or narrower by 
minimizing image energy using an external force and is also influenced by image features such as lines or edges. 
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2.3. Calculation of hippocampus volume on coronal slice 

Volume calculation is done by first calculating the area of each segmented coronal hippocampus 
slice using the active contour method. The real slices of the entire brain MRI image slices are 1-256 slices 
with a thickness, each assumed to be 1 pixel/slice thick. Volume calculation is done by adding up the volume 
of each slice, which is the product of the multiplication of the hippocampus area multiplied by its thickness. 
A total of + 30 slices will be identified for their width, then the volume of each slice is calculated, after which 
the hippocampus volume is calculated from the number of each volume per slice. The volume calculation is 
done using a mathematical approach, namely by adding up each segmented hippocampus slice area, then 
multiplying it by the slice height, which is assumed to be 1 pixel long. Based on the previous calculation, then the 
volume of each slice will be known. Figure 2 shows the plot of the brain MRI hippocampus volume calculation. 
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Figure 1. Coronal slice MRI image 


Segmentation per slice 


Active Contour segmentation 


Calculation ofthe area of each 
slice 


Volume calculation 
[Total number(area per slice x 
height 1 pixel)] 


Figure 2. Flow chart of coronal hippocampus slice volume calculation 


2.4. Three-dimensional visualization 

Biomedical image visualization involves transformation and interaction with multi-dimensional 
biomedical image data sets. Researchers use several methods in creating 3D object views, which can be 
grouped into two techniques: surface rendering and volume rendering [22]. Surface rendering generally 
requires contour extraction, which is the surface of a structure to be visualized and then represented in the 
form of a series of polygons. A tiling algorithm is then applied to place the polygon pieces at each contour 
point to produce a surface appearance. The advantage of this technique is the ability to quickly generate 
image views for a relatively small amount of contour data. Where as the disadvantage is that it lies indiscretion 
when extracting contours and placing polygons on the surface, thereby reducing its accuracy [23]. Volume 
rendering is one of the most powerful techniques in visualization and image manipulation [24]. In this 
technique, there is no need for discretion from the surface so that the integrity of the volume image data is 
well preserved. That is why this technique can produce high-quality displays, although the implication is that 
the computation time is quite long. 
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3. RESULTS AND ANALYSIS 

In this research, we used 40 MRI image data from the OASIS database. The image has identified its 
CDR value. The first step is to cut the image on each MRI image. The cutting was done using free software 
called “MRIcro.” Figure 3 shows the cut image. In coronal slices, one brain MRI consists of + 176 2D image 
slices. The next step is to convert the data file, which initially has the extension “hdr” into a collection of 2D 
coronal brain images with the extension “.bmp.” then the results obtained are as shown in Figure 4, where in the 
image you can see the shape of the brain starting from the first time you can to the end of the 2D image slice. 
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Figure3. Display of MRIcro software 
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Figure 4. Coronal MRI image slice (256 slices) 
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As explained in section 1, this paper will discuss the calculation of the volume of the hippocampus. 
This discussion refers to our previous research [25]. This part of the hippocampus determines the severity of 
Alzheimer’s. Alzheimer’s disease has a small part of the hippocampus, and the size decreases over time. This 
part of the hippocampus was studied starting from the slice when it first appeared to the end of the 
hippocampus slice. In the initial process, we clustered 40 MRI brain image data, which contained a 
hippocampus slice. In the clustering stage, it is done manually, with the eye’s help directly on the parts that 
have the hippocampus in it. We do this because the position of the hippocampus on each slice has different 
coordinates. Also, the hippocampus has a pixel intensity value that is almost the same as the other parts of the 
brain. Figure 5 shows an example of a successfully observed section of the hippocampus. In Figure 5(a) 
shows the temporal horn type of hippocampus, while in Figure 5(b) shows the middle hippocampus and in 
Figure 5(c) shows the hippocampus tail. 

In this experiment, we clustered 40 data on MRI brain images in which there was a slice of the 
hippocampus. We did the clustering process manually by selecting the hippocampus directly. We use the 
manual method because the position of the hippocampus on each slice has different coordinates. It is also 
because the pixel intensity in all areas of the brain has almost the same value. Table 1 shows the results of the 
clusterization in this experiment. 

Referring to Table 1, we see that each brain MRI has a different number of slices. On MRI images 
with a CDR value of 0, the number of hippocampal slices on average is more than CDR 1 and CDR 2. 
Likewise, CDR 1 has more average hippocampal slices compared to CDR 2. The relationship between the 
number of hippocampus slices on the MRI image of Alzheimer’s with the CDR value is like the (1). 


CDR 0 > CDR 1 > CDR2 (1) 


In this experiment, we pre-processed all input data to get a clearer image in distinguishing 
between the hippocampus area and other areas. Here we use the contrast stretching method. In previous 
studies, the value of contrast stretching can be changed in brightness to obtain an optimal coronal 
hippocampus image. However, in this study, we used a controlled variable for the contrast stretching function. 
We do this so that there is no significant difference in pixel intensity values between one image and the next 
image in the initial masking process. Figure 6(a) and Figure 6(b) shows the difference in results between the 
original image and the image after contrast stretching. 

Referring to Figure 6(a) and Figure 6(b), we can see that the image after the contrast stretching process 
has an even distribution of light compared to the original image. The coronal slice of the hippocampus is visible, 
and we can easily carry out the segmentation process. So that with the improvement of image quality using 
contrast stretching, the image becomes more apparent and will facilitate segmentation in the hippocampus area. 

In carrying out the initial masking process, an image enlargement is carried out on the right and left 
hippocampus. This enlargement is done so that the contour points do not hit any part other than the coronal 
side’s hippocampus. Because when the contour points are far from the sides of the hippocampus, the object is 
not appropriately segmented and can affect the volume calculation of the hippocampus. Researchers did not 
use the cropping command in this study because it can change the original Alzheimer’s coronal MRI image 
size. When the size of one slice image is different from the size of another slice image, the 3D visualization 
formation process cannot be processed because 3D image formation requires a 3D matrix of the same size. 
Therefore, the zoom on and zoom off functions are used in the research. 


(a) (b) (c) 


Figure 5. Hippocampus type: (a) hippocampus temporal horn, (b) hippocampus midle, and 
(c) hippocampus tail 
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The active contour process is the following process of initial masking. The segmentation with the active 
contour in this research is based on “edge,” which results in the contours evolving inward so that the contours are 
conical to form the desired object. The working principle of this active contour segmentation is that if the 
external energy from the contour detects a change in intensity (called edge detection), the gradient will reduce 
its speed to close to O so that the internal energy moves along the contour and forms the area of the actual 
object. The result of this active contour segmentation is a binary image, where the hippocampus is 1 (white) and 
other than that it is O (black) with the same size as the original MRI image, which is 208x176 pixels. The results 
of the initial masking will be automatically segmented using the active contour method. Figure 7 describes 
the results of the hippocampus segmentation in this research. 


Table 1. Identification results 


No ID Male/female Age CDR Slicefull Slice hippo 
1  OAS1_0058_MRI1 F 46 0 256 33 
2  OAS1_0044_MRI1 F 47 0 256 30 
3 OAS1_0020_MR1 F 48 0 256 29 
4 OAS1_0011_MRI1 F 52 0 256 31 
5 OAS1_0002_MRI1 F 55 0 256 32 
6  OAS1_0026_MR1 F 58 0 256 33 
7 OAS1_0030_MRI F 65 0 256 31 
8 OAS1_0001_MRI1 F 74 0 256 32 
9  OAS1_0033_MRI F 80 0 256 28 
10 OAS1_0019 MRI F 89 0 256 30 
11 OAS1_0069_MR1 M 33 0 256 25 
12 OAS1_0018_MRI1 M 39 0 256 25 
13 OAS1_0074_MR1 M 43 0 256 31 
14 OAS1_0034_MRI1 M 51 0 256 24 
15 OAS1_0114_MR1 M 62 0 256 34 
16 OAS1_0135_MR1 M 64 0 256 26 
17 OAS1_0130_MR1 M 68 0 256 28 
18 OAS1_0010_MR1 M 74 0 256 32 
19 OAS1_0110_MRI1 M 84 0 256 31 
20 OAS1_0065_MR1 M 90 0 256 30 
21 OAS1_0184_MRI1 F 65 1 256 32 
22 OAS1_0382_MR1 F 67 1 256 27 
23 OAS1_0067_MR1 F 71 1 256 28 
24 OAS1_0056_MR1 F 72 1 256 29 
25 OAS1_0269_MRI1 F 72 1 256 31 
26 OAS1_0316_MR1 F 72 1 256 26 
27 OAS1_0052_MR1 F 78 1 256 26 
28 OAS1_0053_MR1 F 83 1 256 31 
29 OAS1_0035_MR1 F 84 1 256 23 
30 OAS1_0430_MR1 M 71 1 256 29 
31 OAS1_0424 MRI M 75 1 256 23 
32 OAS1_0452_MRI1 M 75 1 256 29 
33 OAS1_0405_MRI1 M 77 1 256 30 
34  OAS1_0268_MR1 M 78 1 256 29 
35 OAS1_0399_ MRI M 78 1 256 29 
36 OAS1_0134_MRI1 M 80 1 256 23 
37 OAS1_0223_MRI1 M 84 1 256 29 
38 OAS1_0031_MRI1 M 88 1 256 27 
39 OAS1_0308_MRI F 78 2 256 22 
40 _OAS1_0351_MRI1 M 86 2 256 23 


Calculation of the volume of the hippocampus can be done when all coronal slice hippocampus 
images have been segmented thoroughly, from the first to the last slice, and entered in the same folder-then 
given the sequential names starting from the beginning of the emergence of the hippocampus (horn) to the 
end of the hippocampus (tail). This sequential naming is used to make it easier to access files in a folder so 
that users do not need to access images one by one, simply by selecting a folder in which there is a collection 
of MRI image slices. The method used to find the volume of the hippocampus is to enter a collection of 
images in the folder into a 3D matrix. The image used is a binary image, so the value of the hippocampus in 
each slice is expressed as 1 (white) and outside the hippocampus is 0 (black). So it can be concluded that the 
volume of the hippocampus is the addition of the 3D matrix into one 1x1 matrix. The 3D visualization process 
is carried out by placing each segmented 2D hippocampus slice into a 3D matrix. In this research, the zero 
function is used to place a 2D image into a 3D matrix. When the 2D coronal Alzheimer’s image is placed on a 
3D matrix, it will look like the illustration visualization in Figure 8. This image is a display of a 2D brain slice 
containing the hippocampus image. Then the images are stacked into a 3D matrix. 


TELKOMNIKA Telecommun Comput El Control, Vol. 21, No. 1, February 2023: 123-132 


TELKOMNIKA Telecommun Comput El Control Oo 129 


In this research, only the hippocampus is shown in 3D, and other parts are ignored. Because the parts 
of the hippocampus in each coronal MRI image have been segmented in the previous processes, only part of the 
coronal hippocampus is automatically displayed on the screen. In the segmented image, the visualization of the 
3D matrix in each stack will resemble the shape of the hippocampus itself, as shown in Figure 9. 

Referring to Figure 9, then we reconstruct the 3D visualization of the previously formed 3D matrix. 
In this research, the results of 3D reconstruction from a collection of 2D MRI images of the hippocampus are 
smoothed using the smooth function; this visual refinement is functional. The image is better than the 
previous results, as shown in Figure 10. Figure 10 shows the results of 3-dimensional visualization, where in 
Figure 10(a) shows the results before smoothing implementation, while in Figure 10(b) is the result of 
visualization after the smoothing process. 


(a) (b) 


Figure 6. Before and after contrast stretching: Figure 7. Segmentation result 
(a) original image and (b) after contrast stretching 
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Figure 8. Alzheimer’s brain coronal slice stacks 


All results of the 3D visualization reconstruction of the coronal slice hippocampus image are shown 
in Table 2. Referring to Table 2, the relationship between volume and CDR shows that the average value for 
CDR 0 on the left hippocampus volume is 1973 pixels, the right hippocampus volume is 2055 pixels, and the 
total hippocampus volume is 4028 pixels. The average value for CDR 1 on the left hippocampus volume is 
1252 pixels, the right hippocampus volume is 1253 pixels, and the total hippocampus volume is 2505 pixels. 
Meanwhile, the average value for CDR 2 on the left hippocampus volume is 731.5 pixels, the right 
hippocampus volume is 671 pixels, and the total hippocampus volume is 1402 pixels. It can be concluded 
that the volume of the left and right hippocampus of each hippocampus is not much different. 
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Figure 9. Hippocampus coronal slice stacks 
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Figure 10. Hippocampus 3D visualization view (a) 3D vizualisation before smoothing and 


(b) 3D vizualisation after smoothing 


Table 2. CDR VS volume of hippocampus 
CDR Left volume Right volume Total volume 


0 1973 2055 4028 
1 1252 1253 2505 
2 731.5 671 1402 


The relationship between CDR and the number of slices resulting from 3D visualization shows that 
the slices that are primarily found in CDR 0 are 30 slices, then CDR 1 is 28 slices, and CDR 2 is 23 slices. 
This shows that the more slices of the hippocampus’s coronal MRI image, the volume will also be more 
significant, as shown in Table 3. In Table 3, we can see that the larger the CDR value, the smaller the number 
of coronal cross-sectional slices produced. 

Meanwhile, the relationship between CDR and user operating time-based on 3D visualization shows 
that at CDR 0, the time used by the user is 74.1 minutes. CDR 1 is 66.7 minutes, and CDR 2 is 53.5 minutes. 
This condition shows that the lower the CDR value of an image, the longer it will calculate the volume. This 
condition happens because the lower the CDR value, the greater the volume and the number of slices, as shown 
in Table 4. 

The relationship between volume and volume of the hippocampus CDR for women and men in our 
3D visualization shows that the average value for the CDR 0 on the female hippocampus volume amounted 
to 3937 pixels, and the male hippocampus volume amounted to 4120 pixels. The average value for CDR 1 on 
the volume of the female hippocampus is 2421 pixels, and the volume of the male hippocampus is 2590 
pixels. Furthermore, the average value for CDR 2 in the volume of the female hippocampus is 1558 pixels, 
and the volume of the male hippocampus is 1247 pixels as shown in Table 5. 

The relationship between volume and CDR from our 3D visualization results for each age shows 
that at CDR 0, the older the Alzheimer’s sufferer, the smaller the volume value. Meanwhile, for CDR 1 
Alzheimer’s sufferers starting from the age of 60 years and over. Meanwhile, CDR 2 for Alzheimer’s patients 
starts at the age of 70 years as shown in Table 6. 
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Table 3. CDR VS the number of slices Table 4. CDR vs time processing 
CDR The number of slces CDR Time processing 

0 30 0 74.1 

1 28 1 66.7 

2 23 2 53.5 

Table 5. Male vs female Table 6. CDR vs age 

CDR Male himppocampus Female hippocampus CDR 40-50 51-60 61-70 71-80 81-90 
0 3937 4120 0 4355 4642 4007 3483 3132 
1 2421 2590 1 - - 3392 2500 2076 
2 1558 1247 2 - - - 1247 1558 


4. CONCLUSION 

3D visualization calculations in the hippocampus coronal slice area have different volume values. 
The average volume value of CDR 0 is more significant than CDR 1 and CDR 2, and the average volume 
value of CDR 1 is more significant than CDR 2 but is smaller than CDR 0. The value of the volume of the 
hippocampus is the sum of the 3D matrix from the 2D image slices of each hippocampus. The results of the 
3D hippocampus visualization on each image also have a different appearance. When segmenting each 2D 
MRI image slice, part of the hippocampus of each image has a different shape and coordinates. The lower the 
CDR value, the bigger the volume and the more slices, and the calculation time by the user will be longer. 
The difference in the volume of the left and right hippocampus did not show a significant difference. 
Regarding gender, the differences in the volume of the hippocampus of men and women are relatively the 
same. As for age, for people with Alzheimer’s, as they get old, the volume will shrink. Moreover, people 
with Alzheimer’s with CDR 1 starting from the age of 60 years and over. Meanwhile, CDR 2 starts from the 
age of 70 years and over. 
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